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We introduce a variation of direct maximal center gauge fixing: the "direct Laplacian" center gauge. The new procedure 
consists of first fixing to the Laplacian adjoint Landau gauge, followed by overrelaxation to the nearby Gribov copy of the 
direct maximal center gauge. Certain shortcomings of maximal center gauge, associated with Gribov copies, are overcome in 
the new gauge, in particular center dominance is recovered. 



1. DRAMA OF GRIBOV COPIES: THE FIRST 
FIVE ACTS 

We proposed a procedure for identifying center 
vortices in lattice configurations based on center pro- 
jection in maximal center (or adjoint Landau) gauge, 
and accumulated evidence in favor of the center vor- 
tex model of color confinement [JxJ] . MCG fixing, 
however, suffers from the Gribov copy problem: The 
iterative gauge fixing procedure converges to a local 
maximum which is slightly different for every gauge 
copy of a given lattice configuration. The problem 
seemed quite innocuous at first; we observed that vor- 
tex locations in random copies of a given configura- 
tion are strongly correlated. Recently, the successes 
of the approach based on MCG fixing have been over- 
shadowed by serious difficulties, so serious that they 
led the authors of Ref. [[|] to speak about the Gribov 
copy drama. Indeed, it follows the structure of the 
classical ancient tragedy: 

1 . Exposition: The projected string tension in MCG 
was shown to reproduce the full asymptotic string ten- 
sion (center dominance) [jj]]. 

2. Complication: Kovacs and Tomboulis ob- 
served that if we fix to Landau gauge, before over- 
relaxation to MCG, center dominance is lost! 

3. Climax: Bornyakov et al. [^] pointed out strong 
dependence of results on the number of gauge copies 
used in MCG maximization. They chose the best 
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from N CO p copies and extrapolated for N cop — > °°: the 
full string tension was underestimated by about 30% ! 

4. Reversal: It was shown ^\ that the results de- 
pend on the number of gauge copies and the lattice 
size. On a large lattice (compared to the typical size 
of the vortex core) the problem disappears. 

5. Catastrophe: Bornyakov et al. observed that 
using simulated annealing one can find better maxima 
than with the usual method of overrelaxation, but the 
center projected string tension again is only about 2/3 
of the full string tension. 

Our aim in this paper is to argue that, unlike an- 
cient tragedies, the Gribov copy drama can have the 
sixth act, a happy ending. We shall introduce a new 
procedure, direct Laplacian center gauge (DLCG), 
that overcomes the difficulties associated with Gribov 
copies. Before doing that, however, we shall briefly 
review a new insight into MCG fixing, due to 

2. GAUGE FIXING AS FINDING A "BEST FIT" 

Imagine running a Monte Carlo simulation and ask- 
ing for the pure gauge configuration closest, in con- 
figuration space, to the given thermalized lattice. It is 
easy to show that finding such an "optimal" configu- 
ration is equivalent to fixing to the Landau gauge. 

Allow now for Zi dislocations in the gauge trans- 
formation, i.e. fit the lattice configuration by one with 
thin center vortices: 

U™{x) = g(x)Z,(x)g\x + jj), Z fJ (x) = ±l. (1) 
Such U™ r (x) is a continuous pure gauge in the ad- 
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joint representation, which is blind to the ZJx) factor. 
Then the fit can be easily performed in two steps: 

1. one can determine g(x) up to a Z2 transforma- 
tion by minimizing the square distance between Ua^ 
and U A ™~ in the adjoint representation - this is equiv- 
alent to fixing to the adjoint Landau gauge, aka direct 
maximal center gauge; 

2. Zu(x) is then fixed by the center projection pre- 
scription: Z M (x) = sign Tr [g^ (x)U fl (x)g(x +/>)] . 

It is now immediately clear that UJ l " r (x) is a bad 
fit to UJx) at links belonging to thin vortices (i.e. to 
the P-plaquettes formed from Z jJ (x)). We recall that 
a plaquette p is a P-plaquette iff Z(p) = — 1 (where 
Z(C) denotes the product of ZJx) around the con- 
tour C) and that P-plaquettes belong to P-vortices. Let 
us write the gauge transformed configuration as 

g U M (x) = Z^ (x) e iA »^ , Tr e iA " w > 0. (2) 

At large (3 values jTrfl/p] = 1 — O(jj), and equals to 

(Zp) ±TrV\e iA ^ x) = (-nxiTrTTe^- (3) 
A f A onP-plaq. z A /) A 

The last equation implies that at least at one link 
belonging to the P-plaquette Af,(x) cannot be small, 
therefore g U^(x) must strongly deviate from the cen- 
ter element. This indicates that the quest for the global 
maximum may not be the best strategy; one should 
rather try to exclude contributions from P-plaquettes 
where the fit is inevitably bad [J7j], or modify the gauge 
fixing procedure to soften the fit at vortex cores. 

3. DIRECT LAPLACIAN CENTER GAUGE 

Our new proposal to overcome the Gribov copy 
problem was inspired by the Laplacian Landau gau- 
ge The idea is the following (for details see [|[|): 

To find the "best fit" to a lattice configuration by a 
thin center vortex configuration one looks for a matrix 
M{x) maximizing the expression: 

%m = E Tr \ mT {x)UAn{x)M{x+ii)\ , (4) 

with a constraint that M(x) should be an SO(3) matrix 
in any site x. We soften the orthogonality constraint 
by demanding it only "on average": (M T ■ M ) = 
±Z.M T (x)-M(x) = l. 

It is convenient to write the columns of M(x) as a 
set of 3-vectors f a {x): fa{x) —M a j,(x). The optimal 



M(x) is determined by three lowest eigenvectors fj(x) 
of the covariant adjoint Laplacian operator: 

2>y (*,y) = 2D8, T 5 i7 - £ [U A . ±M (x)} 5 y , x±[l . (5) 

The resulting real matrix field M(x) has further to 
be mapped onto an SO(3)-valued field gA(x). A naive 
map (which could also be called Laplacian adjoint 
Landau gauge) amounts to choosing gA (x) closest to 
M(x). Such a map is well known in matrix theory and 
is called polar decomposition. 

A better procedure, in our opinion, is the Laplacian 
map, that leads to direct Laplacian center gauge. We 
try to locate gA(x) as close to M(x) local maximum 
of the MCG (constrained) maximization problem. To 
achieve this, we first make the naive map (polar de- 
composition), then use the usual quenched maximiza- 
tion (overrelaxation) to relax to the nearest (or at least 
nearby) maximum of the MCG fixing condition. 

4. CENTER DOMINANCE 

To test the procedure of the last section, we have 
recalculated the vortex observables introduced in our 
previous work with P-vortices located via cen- 
ter projection after fixing the lattice to the new direct 
Laplacian center gauge. The full set of results has 
been published in [gp. Here we just concentrate on 
the issue of center dominance which is crucial for the 
whole picture. Without being able to reproduce the 
string tension of the full theory, one cannot claim to 
have isolated degrees of freedom that are associated 
with the confinement mechanism. 

A subset of our results, obtained at a variety of cou- 
plings on our largest lattice sizes, is shown in Fig. [j]. 
It is clearly seen that center dominance is restored 
in DLCG Moreover, we once again observe preco- 
cious linearity (the very weak dependence of pro- 
jected Creutz ratios on the distance R) signalling that 
the center projected degrees of freedom have isolated 
the long range physics, and are not mixed up with ul- 
traviolet fluctuations. Creutz ratios differ by at most 
about 10% from the phenomenological value of the 
string tension. 

Also other physical results are reproduced in the 
new gauge: scaling of the vortex density, agreement 
of ratios of vortex limited Wilson loops with simple 
expectations, the fact that removing center vortices 
from lattice configurations destroys confinement [H]. 
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DLCG, Center Dominance and Precocious Linearity 



Creutz Ratios at p=2.3 from Products of Z 2 Loops 
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Figure 1. The ratio of projected Creutz ratios to the 
full asymptotic string tension, as a function of R. 
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Figure 2. Creutz ratios from "product Wilson loops' 
show strong correlation between DLCG and ILCG. 



5. COMPARISON TO INDIRECT LCG 

To avoid problems with Gribov copies, de Forcrand 
et al. [ |lO| ] proposed the Laplacian center gauge. The 
procedure is similar to ours, but differs in a few im- 
portant points. The gauge transformation matrix is 
found from two lowest eigenvectors of the covariant 
adjoint Laplacian operator, one first fixes SU(2) to the 
U(l) subgroup, and center vortex surfaces should in 
principle be identified via ambiguities of the gauge 
fixing condition. However, there is no good sepa- 
ration between confinement and short range physics 
in this gauge (center dominance seen only at largest 
distances, no precocious linearity, vortex density not 
scaling), and identification of vortices on a lattice 
via gauge fixing ambiguities is practically impossible 
(center projection is necessary). Langfeld et al. flTT| ] 
therefore used overrelaxation to MCG after fixing to 
LCG. This we would caU indirect LCG (ILCG). 

It turns out that results from DLCG and ILCG are 
quite similar. Center dominance, precocious linear- 
ity, and vortex density scaling are observed, how- 
ever, Creutz ratios in ILCG are somewhat lower than 
those in DLCG. One can attribute the similarity of 
both procedures to the strong correlation between P- 
vortex locations in both gauges. This is illustrated in 
Fig. ^| showing Creutz ratios calculated from "prod- 
uct Wilson loops" W prod (C) = (Z IWG (C) Z DWG (C)). 
If there were no correlation, these ratios should ap- 
proach twice the asymptotic string tension, while they 
should go to zero for strong correlation. Our data 
clearly favor the latter case. 



6. CONCLUSION 

We have proposed a new gauge, DLCG, that com- 
bines fixing to Laplacian adjoint Landau gauge with 
the usual quenched maximization. The first step of 
the procedure is unique, in the second step no strong 
gauge copy dependence appears. The procedure can 
be interpreted as a "best fit" of a lattice configura- 
tion by thin vortices, softened at vortex cores. Cen- 
ter dominance, precocious linearity and scaling of the 
vortex density are recovered in the new gauge. 
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